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Abstract 

We show that Einstein's equations in a nonstandard physical gauge 
have vacuum solutions with an asymptotically fiat rotation curve as it 
is observed in the dark halos of galaxies. Introducing a material disk 
into this model we find a matter density in accordance with the Tully- 
Fisher relation. 
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1 Introduction 



General relativity is a classical gauge theory. This implies that the funda- 
mental fields, like the metric g a a, are not directly observable. Therefore, in 
investigating gravitational effects it is important to identify the observable 
quantities which are actually measurable. Every observable is defined by a 
measuring process. On the scale of galaxies one most important observable 
is the circular velocity V(r) of stars or gas which can be measured by the 
Doppler shift of spectral lines (r is the radius of the circular orbit). 

In the theoretical analysis one should try to relate the observables to the 
fundamental fields, say to g a p{x). The best would be if the metric under 
certain assumptions can be uniquely expressed by the observables. Then the 
gauge ambiguity has been removed, the gauge is fixed by a physical require- 
ment. The other possibility is to choose the gauge on unphysical grounds, 
for example by some geometric convention or / and to simplify the solution of 
the differential equations. This is dangerous because one might miss some 
important physics. Our program of fixing the gauge by observables is a sort 
of inverse procedure compared with standard general relativity where one 
first calculates a metric by solving Einstein's equations in some special gauge 
and then determines the observables. Clearly, in standard general relativity 
one cannot be sure that one finds all physically relevant solutions. Indeed 
we are going to show that Einstein's equations have vacuum solutions with 
an asymptotically flat rotation curve V(r). These nonstandard solutions 
can be used to describe the dark halo of galaxies without introducing hy- 
pothetical dark matter. After the recent Xenon-experiment has again not 
found any signal of dark matter particles [1] one should seriously investigate 
nonstandard general relativity. 

The paper is organized as follows. As preliminaries we first consider 
the rotation curve in a general spherically symmetric metric. Although this 
might be widely known we do not know a good reference. In Sect. 3 we solve 
the vacuum Einstein's equations in our general spherically symmetric setting 
and express the metric tensor by the circular velocity V(r). From these 
vacuum solutions we construct in Sect. 4 solutions with a disk of ordinary 
matter by means of the well-known displace, cut, and reflect method. This 
gives a simple model of a spiral galaxy. Assuming a circular velocity which 
is constant = Vflat for large r, we find a matter density proportional to Vfl at . 
This is in accordance with the Tully-Fisher relation [2] [3]. 
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2 The circular velocity in general relativity 



We consider a star moving in an arbitrary static gravitational field. Follow- 
ing Weinberg ([4], chapter 3/2) we introduce the freely falling coordinate 
system £ Q of the moving star. In this system the star is at rest and the 
observer on earth moves with the 4-velocity 

u a = ^f- = (u°,H), u° = (1 - V 2 )- 1 / 2 , u = u°V, (2.1) 

CLT 

because in the locally inertial coordinates special relativity holds. Here r is 
the proper time 

dr 2 = Va ^ a d^, (2.2) 

Va/3 = diag(l, —1,-1,-1) is the Minkowski tensor and the speed of light is 
put = 1. In addition Weinberg introduces a laboratory coordinate system 
x M which in our case is attached to the observers telescope. To simplify the 
following discussion we assume that the astronomer on earth has corrected 
his measurements for the motion of the earth with respect to the center 
of the galaxy, so that we can choose the center of the galaxy as origin of 
the laboratory coordinate system. Now the star moves on a geodesic ([4], 
eq.(3.2.3)) 

= » (2-3) 

where T^o are the Christoffel symbols of the metric g^ v . The latter is defined 
by ([4], eq.(3.2.7)) 

Since the local inertial coordinates £ a can be changed by arbitrary Lorentz 
transformations, we can choose a Lorentz boost such that the two coordinate 
systems are at rest with respect to each other. Then we have 

£ = 0, and ^ = (2.5) 
for j = 1, 2, 3. This leads to 

500 = teJ ' 9 * = -±<WdxS (2 - 6) 

1=1 
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and goj = 0. Now it follows 

o de dfdx dx° 

u= ^ = WlF = ^lt- (2 - 7) 

In the following we assume that the non-diagonal elements gjk vanish. 
Then from the invariant 

B , Ch2 ^2 n di a d^dx^dx" 

u W = ^ 2 - u2 = 1 = B£*#top-&F-&F = 

dx 11 dx v ( dx° \ 2 / dxi \ 2 

we find 

f /2 __J_ n f dxi \ 2 _ 9jj (dxi/dr) 2 

V ~ (u^ 9j A dT J- ^ 500 (dx°/dr) 2 ^ 

where (2.7) is used. 

We want to specialize this for circular motion r = const. Now the 
geodesic equation (2.3) for /i = 1 reads 

Taking the circular orbit r = const in the equatorial plane # = 7r/2, then 
(2.9) gets simplified to 

and from (2.8) we finally obtain the important relation 

ffoo r 33 

In the following we consider static spherically symmetric metrics of the 
form 

ds 2 = g^dx^dx" = e a dt 2 - e b dr 2 - rV(cM 2 + sin 2 tf# 2 ), (2.12) 
where a(r), b(r), c(r) are functions of r only. For this metric we have [5] 



r 1 - —p a - b 
1 oo — - e 
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rL = -(yc' + r)e c - b 
here the prime means d/dr. Then (2.11) simply becomes 



(2.13) 




d + 2/r' 



(2.14) 



3 Solution of the vacuum equation 

The metric functions a, b, c appearing in (2.12) must satisfy differential equa- 
tion which follow from Einstein's equations. In standard general relativity 
one puts c = 0. This is a special choice of gauge which leads to Birkhoff's 
theorem and the Schwarzschild metric. This works well in the solar system, 
but obviously not on the galactic scale. The standard way out is to abandon 
the vacuum equations and assume some hypothetical dark matter. As long 
as this dark matter is not convincingly recorded one should study the other 
possibility of retaining c(r) 7^ 0. Then the vacuum solution is no longer 
unique. To fix it uniquely we take the expression (2.14) for the circular 
velocity V(r) as our nonstandard gauge condition. In other words, every 
circular velocity V(r) defines a particular gauge. That means V{r) must be 
given, so that the theory seems to have less predictive power. What seems 
to be a weakness is a strength: The asymptotic V(r) cannot be predicted on 
the basis of the vacuum equations alone, the dynamics of the normal matter, 
that means the detailed structure of the galaxy, must necessarily be taken 
into account. Indeed a uniform asymptotic velocity profile for all galaxies 
seems not to exist. Only with c 7^ is it possible to carry out our program 
to express the metric by the observable V(r). 

The non- vanishing components of the Ricci tensor for the metric (2.12) 
are the diagonal elements [5] 



Rtt 




a'b' + a'c' + -a') 



(3.1) 



r 
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the prime always denotes d/dr. Then the Einstein's equations without mat- 
ter can be reduced to the following three differential equations 



G f , = e 



a—b 



~ c " ~ 7 C ' 2 + l b ' c ' + V " 3c ')l +\(e a 
4 2 r J r 2 



c - e a - b ) 



(3.5) 

IaV + V + c/) + ^ + l(l-c*- c )=0 (3.6) 



Gaa — — e 



r 2 „_ b 



a » +c »A( b >- a '-2c')+-(a' 2 -a'b'+a'c'-b'c'+c' 2 ) = 0. (3.7) 
r 2 



It is not hard to see that there are only two independent field equations. 
Indeed, using (3.6) b can be expressed by a and c. Eliminating b in (3.5) 
and (3.7) there results one second order differential equation for a and c: 



n a " ( i 2 \ 4 ii 2. , 

c" = — c' + - +-7 + 0V + - + -(a' + c' . 
a / \ r j t a 2 r 

Introducing the new metric function 

f(r) = c(r) + 2 log - 



c' 2 



(3.8) 



(3.9) 



where r c has been included for dimensional reasons, equation (3.8) assumes 
the simple form 

f" a" f 

T-H = J + T (3 ' 10) 



This can immediately by integrated 

/' / 

log — = a H h const. 

a' 2 

On the other hand the circular velocity squared (2.14) becomes 

V 2 (r)^u=f r 



(3.11) 



(3.12) 



It is the velocity squared u(r) which appears in all equations. Using (3.12) 
in (3.11) we have 

a+ £ +log— = -log it. (3.13) 

2 r c 

Differentiating and eliminating d by means of (3.9) and (3.12) in the form 

2 a' 2 



c' = r- 



u 



(3.14) 
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we find 

In 1 

= ^_ (3 15) 

1 + 2u V ' 

This gives the first diagonal element of the metric 

gtt = e a = — ^- (3.16) 
m l + 2u v ' 

where K a is a constant of integration. Then from (3.13) we get 

5^ = e c = ( ) 3.17 

where -ftT c is another integration constant which contains r c . Finally, expft 
follows from (3.6) 

9rr = e b = K c (^f{l + 2u). (3.18) 

We have succeeded in expressing the metric by the circular velocity squared 
u(r). If we choose u = r s /2(r — r s ), we recover the standard gauge c = 
and the Schwarzschild metric. But now also flat rotation curves are possible 
without dark matter. As the ether is superfluous in special relativity, so 
seems to be dark matter in general relativity. From (3.16-18) we are able to 
predict other observable quantities which can be computed from the metric, 
for example lensing data [5]. In this way the theory can be tested. Another 
test is investigated in the next section. 



4 Thin material disk with a dark halo 

We study a simple model of a spiral galaxy by assuming that the normal 
matter is concentrated in the equatorial plane z = with a singular density 
~ 5 1 (z). For this problem the theory of distribution valued curvature tensor 
is appropriate which is mainly due to Taub [7] . To be self-contained we give 
a simple derivation of the relations we need. Let 5 be a three-dimensional 
surface in 4-space where the metric tensor g„ v is continuous but has finite 
jumps in the normal derivatives; the derivatives in the tangential directions 
are assumed to be continuous. In an admissible coordinate system let S be 
described by the equation 

<p(x) = (4.1) 



and have the normal vector 
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Then the finite discontinuities in the first partial derivatives of g^ u are given 
by 

[9 ^°i^\ + - d i^\= n ^ (4 - 3) 

where + and — mean the limiting values from both sides of S. This follows 
from the decomposition of the gradient into normal and tangential compo- 
nents. The corresponding jumps in the Christoffel symbols then are 

2[T%\=npb« + n 1 b<$-n a bp 1 . (4.4) 

The Ricci tensor 

R/xv = d a r° v — d u r^ a + r^r^ — r^r^ (4.5) 

contains derivatives of T, consequently the finite jumps lead to singular 
contributions proportional to the delta distribution 8s with support on S 
according to 

9^1^ = F^npSs. (4.6) 
Then it follows from (4.4) that 

^Vising = ^("""""V + n< *V*nz/ + n a b au n^)5s, (4.7) 

with 

y*p = b a p-\b8%, b = <rb liV . (4.8) 

This is in agreement with eq.(2.14) of Taub, note that his convention for the 
Ricci tensor is the negative of our (4.5). 

In the Einstein's equations these singular distribution must be compen- 
sated by a distribution valued energy-momentum tensor 

(Rfiv - Rising = Kt^vSs, (4.9) 

where 

R = g^R a p, « = ^- (4.10) 

If the jumps b^ u of the normal derivatives of g^ v are known, can be 
calculated from (4.7) and (4.9): 



o . 2 ( i a n n ii\( t n n v 

-2«V = n ^ - —f-)(g v ~ —J-) 
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where n 2 = n a n a . This agrees with eq.(6-2) of Taub. The singular contri- 
bution (4.11) must be added to the regular energy-momentum tensor which 
renders the field equations fulfilled outside of the surface S. 

Now we come to our simple galaxy model where the normal matter is 
concentrated in the plane z = which is our singular surface S. Outside this 
plane we have vacuum with a dark halo as it is described by the nonstandard 
spherically symmetric solution (3.16-18). To have a simple representation 
of the plane z = and the corresponding delta-measure we go over to 
cylindrical coordinates (t, R, z, (f) 

R 

r 2 = R 2 + z 2 , z = rcos'd, sini? = — . (4.12) 

r 

Then the metric (2.12) assumes the following non-diagonal form 

ds 2 = g^dx^dx" 

with 

a R b Z c c Z b 

5oo = e , 5n = T e ~e , g 2 2 = ~-e ~e 

j^>Zi ry* £j ry* A y 

gn = 921 = -2^(e b - e c ), g 33 = -R 2 e c . (4.13) 

For simplicity we still write r, but our admissible coordinates are x 1 = 
R,x 2 = z. We also need the inverse 

9 m = e~\ 9 U = 9 -§, 9 22 = f 

9 12 = - 9 -§=9 2 \ 9 33 = ~, (4-14) 
L> 533 



where the determinant D is equal to 



D mifl-n - (fin) 2 --< Jl+l -3^(e b -e c ) 2 . (4.15) 



To construct the metric with the material disk we apply the widely used 
displace, cut, and reflect method which goes back to Kuzmin [7] and since 
then has been used by many authors. Following the procedure of Voigt and 
Letelier [8] we take the metric (3.7) in the half space z > d > 0, displace it 
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to z = and reflect it for z < 0. This produces the finite jumps in the z- 
derivatives of g^ v . The whole procedure is equivalent to the transformation 
z — > \z\ + d. The normal vector is = (0, 0, 1, 0) = 5 2 and 

v vu v2 v 22 

n = g ^ = g , n n v = g . 

The jumps (4.3) in the normal derivatives on z = which we need are equal 
to 

611 = [5ii,2] =9[™-f 2 e c (4.16) 

&33 = [533,2] =933 — ' 

where the prime always means d/dr keeping z and R constant. Now from 
(4.11) we find the energy density 

* = sK + ^ i ») (4 ' 17) 

with D = 511522 - 9i2- Using 

2dr , 2 c x 2d , 

&11 = — ( 5n el, 633 = — 533c , 

f \ r ' t 

we finally obtain 

4 = --(^^(4) + %^ 2c - ^ e ~ c )- ( 4 - 18 ) 

nr \ r 4 r z r s R A > 

Here we have to put z = everywhere which gives r 2 = R 2 + d 2 . 

Now we must specify the circular velocity squared u(r) in order to fix 
the metric. We are particularly interested in the case of an asymptotically 
flat circular velocity which in the usual terminology corresponds to a dark 
halo. Therefore we assume u(r) of the form 

n(r) = Uflat + ^ + 0(r- 2 ) (4.19) 
r 

for large r. Then it follows from (4.16-18) 

e a = K a + 0{r-\ e b = ^+0(r- 5 ) 



e c = ^ + 0(r- 3 ) (4.20) 
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(4.21) 




2d r 2 
~kL c ~R? 



(l + 0(R- 1 )). 



(4.22) 



This is proportional to the density of normal matter because we consider a 
static energy-momentum tensor. Taking (4.21) into account we find that 



for large R. This is in accordance with the baryonic Tully-Fisher relation for 
galaxies [2] [3], which states that the total baryonic mass M is proportional 
to t . In fact, one can show [10] that the contribution of the inner part 
R < Ri of the disk can be made arbitrarily small compared to the outer 
part between R\ < R < R2, say. We emphasize that M is obtained from i§ 
by integrating with the Euclidean surface measure R dR dcf), because this is 
what astronomers are doing when they determine M from luminosity mea- 
surements. Our theory gives a very natural explanation of the Tully-Fisher 
relation which , otherwise, theoretically and observationally is somewhat 
mysterious. 

The radial pressure t r r vanishes because G rr (3.6) does not contain a 
second derivative. Therefore our model must be interpreted as a dust disk 
with purely azimuthal stresses. This is not very realistic and it remains to 
be investigated whether the Tully-Fisher relation is a generic property for 
more physical galaxy models. 
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